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Abstract 



The magnetic dipole (Ml) and electric quadupole (E2) responses of two- 
dimensional quantum dots with an elliptic shape are theoretically investigated 
as a function of the dot deformation and applied static magnetic field. Ne- 
glecting the electron-electron interaction we obtain analytical results which 
indicate the existence of four characteristic modes, with different S-dispersion 
of their energies and associated strengths. Interaction effects are numerically 
studied within the time-dependent local-spin-density theory, assessing the va- 
lidity of the non-interacting picture. 

PACS 73.20.Dx, 72.15.Rn 
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I. INTRODUCTION 



The electric dipole (El) response of quantum dots has deserved much attention in recent 
years, mainly motivated by the measurements of far-infrared absorption in these systems 
IHII. The manifestation of the so called magnetoplasmons in a perpendicularly applied mag- 
netic field constitutes an example of collective oscillation in finite Fermi systems, physically 
analogous to those existing in metal clusters, atoms or nuclei. In parabolically confined dots 
magnetoplasmons are understood as rigid motions of the electronic center of mass, a result 
that stems from the generalized Kohn's theorem However, deviations from parabol- 

icity, such as angular deformations or non-quadratic radial behaviour, may result in more 
complicated absorption patterns which are currently being much investigated 



The understanding of other excitation modes of quantum dots are also essential for 
a proper characterization of these systems. Techniques based on resonant inelastic light 
scattering have already proved extremely useful to this purpose |IT|-p^[. In fact, using 



polarization selection rules they permit to disentangle charge- density from spin-density and 
single-particle excitations, as well as to discern different multipolarity peaks in each channel. 
The momentum transfer and magnetic field dependence of the different excitations of a single 
sample have been recently reported using this technique WM. It is worth to mention that 



theoretical analysis of the Raman spectra in circularly symmetric quantum dots have been 
presented in Refs. [|T^-[I7|]. 

For elliptical quantum dots orbital-current modes and their relationship with the 
quadrupole response were analysed by us in Ref. [ITSJ. A characteristic low energy mode, 



depending on deformation and with a conspicuous signal in the magnetic dipole (Ml) chan- 
nel was found in Ref. [113, while the relevance of orbital excitations in the build-up of the 



electronic moment of inertia was discussed in Ref. ||19|| . In this paper we focus on the mag- 
netic field dependence of the orbital and quadrupole modes, a necessary step to understand 
the Ml and E2 responses in deformed dots and towards the description of Raman scatter- 
ing in symmetry unrestricted nanostructures. It will be shown that four excitation modes 
with 5-dependent energies and strengths characterize both Ml and E2 reponses in ellip- 
tic dots. This result follows from analytical calculations in the so-called non-interacting 
deformed-harmonic-oscillator (NIDHO) model, as well as for interacting electrons in the 
local-spin-density approximation (LSDA). Therefore, it constitutes a robust picture against 
interactions in a somewhat similar way as Kohn's modes are. 

This paper is organized as follows: Sec. II is devoted to the NIDHO model for Ml and 
E2 responses; in Sec. Ill we present the time-dependent LSDA results for the same channels; 
finally, Sec. IV presents the conclusions. 

II. THE NIDHO MODEL 

A. The analytic solution of the ground state 

We consider a system of non-interacting electrons, confined to the xy plane by an 
anisotropic parabola 



\x,y) = l-(u J y + ujy). (i) 
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The mean parabola coefficient ojq = (uj x + u> y )/2 is parametrized as usual, in terms of the 
Wigner-Seitz radius r s and electron number N as Uq = In the following we shall 

label the dot deformation with the parameter 5 = uj v /uj x (ranging from to 1). Taking into 
account a uniform magnetic field B = Be z , the system Hamiltonian reads 

N n i 

H = E ^ (Pi + 7A(r,)) 2 + v^fe) + g* 7 BS z , (2) 
i=i 11 J 

where 7 = e/c (we assume Gaussian magnetic fields) and within the symmetric gauge 
A(r) = B/2(—y, x). The last piece is the Zeeman term, depending on the total spin S z and 
the effective gyromagnetic factor g* ||21|| . Obviously this is a one-electron picture, in which 
the relevant one-electron Hamiltonian is h = h xy + g* / yBs z . Because of the magnetic field, 
the spatial part h xy deviates from a simple harmonic oscillator problem. Namely 



h xy = \{vl + P 2 y ) + ^(xPy - VPx) 
1 



2 v x ry, 2 



+ ^y+uy), (3) 

with uj^. — u^.+ and ujy = Uy + the new parabola coefficients, given in terms of the 
cyclotron frequency oj c = eB/c. 

The nontrivial problem posed by Eq. (|3|) has been elegantly solved in an analytical way, 
by Dippel et al. |22j and Madhav and Chakraborty |23|] in the context of atomic physics and 
quantum dots, respectively. In the following, we will refer to the derivation by Dippel et 
al.. They introduce a similarity transformation to a new Hamiltonian h% = U~ x hX) ', where 
jj _ e iaxy e tpp x py _ parameters a and (3 are chosen in order to obtain two decoupled 

oscillators in the x and y directions 

p 2 x p 2 y M x 2 2 M 2 2 2 



2M a 2M 2 2 1 2 



The reader is addressed to Ref. |22| for the detailed expressions of a, (3, Mi, M 2 , co± and u> 2 
in terms of the original Hamiltonian parameters. 

Most importantly, we can now express the eigenstates of h xy in terms of those of h^. The 
latter are simply products of one-dimensional oscillator functions with a given number of 
quanta, i.e., 

nin 2 

®mn 2 ( x , V) = 0ni (x)<Pn 2 (y) , (5) 

while the eigenstates of h xy are given by 

hxy^f nmal^ J V) B nin2 n\U2 \ x t y) ) 

*nm 2 (x, y) = U $ nin2 (x, y) . (6) 
Notice that the energy eigenvalues of h xy coincide with those of h 3 and that they are 

E nin2 = (ni + + (n 2 + -)o; 2 . (7) 

Having the one-electron states it is then a simple matter to obtain the N electron ground 
state of the NIDHO model just by filling the N lower energy orbitals. 
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B. Ml response 

The excitations within the NIDHO model will correspond to independent particle-hole 
(ph) transitions in the level scheme of the ground state. Taking the orbital angular momen- 
tum l z as excitation operator the Ml strength function will be given by 

Sui(u) = £//»(!-/„) |(nMl^ _1 4C/|n?nf}| 2 x 

hp 

5(E P -E h -u), (8) 

where we have denoted the electron states $ nin2 as |n 1 n 2 ) and the occupation numbers are 
given by the / factors. 

To evaluate the matrix element in Eq. @ we first need to obtain the transformed opera- 
tor. This can be accomplished by using the Baker-Haussdorf lemma for unitary transforma- 
tions [p4" |, which introduces a nested conmutators expansion. A straightforward calculation 
yields 

U-\U = (1 - 2a(3)£ z + (3(1 - a(3)(p 2 x - p 2 y ) + a(x 2 - y 2 ) . (9) 
Position and momentum operators may now be transformed to creation and anihilation ones 

x = ( ^ (at + ai) 

Px = — W - ai) , (10) 



with similar relations for y and p y in terms of a 2 and a 2 . The matrix elements are therefore 
written in terms of the creation and anihilation operator ones 



(n a \a + \n p ) = ^Jnp + 18 nanf}+1 
(n a \a\np) = ^Jn^K a n p -\ ■ (H) 

From the calculation it emerges that there are only four allowed transition energies: 
u)\x = 2ui, u 22 = 2u2, uj + = ui + 0J2 and = \u\ — uj 2 \ whose Ml strengths are 

q ( \ (<x-0(l-<x/3)Mfaf) 2 

SMl(Wil) = 4M 2 U 2 ^ (l-/ni+2n 2 )(^l + l)K+2) 

1 1 occ. levels 

a , \ (a- (3(1- a(3)M 2 ul) 2 

SmiKO = S 4mC 2 ^ ( fnm 2+ 2)(n 2 + l)(n 2 + 2) 

2 ^2 occ . levels 

J2 I 1 - /m+ln 2 +l)(^l + l)(^2 + 1) 
occ. levels 

„ , x (1 -2a(3) 2 , %r w x2 

£ (l-/», +I » a -i)("i + l)n2. (12) 

occ. levels 
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Result fll2|) provides a clear interpretation of the different modes: Wu corresponds to the 
absorption of two x-oscillator quanta, 0^22 absorbs two y quanta while ui + takes one quantum 
in each oscillator. By contrast, u_ is associated with the absorption of one x quantum 
accompanied with the emission of a y quantum. We remind the reader that our convention 
(uj y < uj x ) renders the absorption of a y quantum accompanied by the emission of an x 
quantum energetically forbidden. 

Figure 1 depicts the absorption energies and associated intensities ujS(uj) as a function 
of the magnetic field for different deformations in a quantum dot with N = 6 electrons and 
r s = 1.51. For 5 > 0.7, when the system is just slightly deformed (5 = 1 corresponds to the 
circular case), only the u + mode is active at low B's. As the magnetic field is increased, very 
rapidly the un and U22 modes gain strength and disperse in energy. Eventually, for very 
large magnetic fields U22 becomes the lowest energy mode and carries most of the strength. 

The u;_ mode deserves a special discussion. Deformed dots at low deformation and mag- 
netic field are characterized by the same oscillator occupancy numbers of the corresponding 
circular limit. In the case of magic number dots (N = f\f(J\f + 1) with M an integer number 
specifying the last occuppied shell) the An\ = 1, Ari2 = — 1 ph transitions are Pauli-blocked 
thus completely inhibiting the uj- mode. This excitation is switched on (with a sudden 
increase in strength) only when the deformation or B are high enough to break the closed 
shell structure, which takes place in first instance when 



Figure 2 displays the regions in the (S, B) plane in which the u>- mode is active, for different 
dot sizes. This effect is absent for non-magic electron numbers since there is no Pauli blocking 
for these systems and, as a consequence, U- is active at all 5 and B. For N = 12, 20, . . . , 
subsequent changes in the ground state level structure give also rise to enhancements of the 
strength by allowing an increasing number of ph transitions. The mechanism discussed 
here is also responsible for strong enhancements of the moment of inertia in these systems 



The general situation in Fig. 1 is qualitatively very similar for the deformations 5 = 0.9 
and 0.7, we remark however that for 5 = 0.7 the splitting uon — U22 is higher at low B and 
the cu_ mode switches on at a smaller magnetic field. For large deformations (5 = 0.4) the 
uj_ mode is active even at B = 0. Notice also that the splitting a>n — uj 2 2 is so high that co> 2 2 
approaches uo_ at vanishing magnetic fields. A common feature to all deformations is that 
U22 dominates the response at high B's. Altogether, the NIDHO model provides a quite 
interesting scenario of the Ml channel, with the emergence of novel modes and big strength 
transfers between them as a function of the magnetic field and deformation. 



El- 



C. E2 response 



The E2 strength is given in terms of the quadrupole operator xy by 



Se 2 M = £A(1-/ p ) (n'ln^\U- 1 xyU\n p 1 n p 2 ) x 



hp 




(14) 
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In a way similar to the Ml analysis we may obtain the E2 matrix elements. After a straight- 
forward calculation one finds the same four modes of the Ml channel manisfesting in the 
E2 spectra, although with different strengths (see third column in Fig 1). The detailed 
expressions are 

P 2 

5e2Ki) = - ( 1 -/n 1 +2n 2 )(^l + l)(^l+2) 

occ. levels 

SebM = -J J2 (l-/nin 2 +2)h + l)h + 2) 
occ. levels 

_ , . _ (l-J 2 M 1 M 2 u 1 u 2 f 
Se2{uJ+) ~ M X M^ 2 X 

J2 (l-/n 1+ in 2 +i)(ni + l)(n 2 + l) 

occ. levels 

_ , > (1 + /? 2 M!M 2 ^ 2 ) 2 

^ > E2\^—) — X 

M1M2U1U2 

Yl i 1 - /ni+ln 2 -l)(«l + 1>2 • (15) 
occ. levels 

We notice that the share of strength as a function of the magnetic field and deformation is 
similar to the magnetic dipole case. Several differences may, however, be remarked. Firstly, 
the u> + mode is relatively enhanced for all 5 and B as compared to the results for the Ml 
channel. For instance the intensity ratio 7e2(^+)/^E2(^-) is enhanced by a factor « 4 at 
5 = 0.4 and B = 0. It is also worth to point out that the U\\ and intensities are reverted 
with respect to the Ml result. Finally, at large magnetic fields the low energy 1022 mode is 
no longer the dominant one since both lo\\ and uj- show a higher intensity. 

As remarked above the number of electrons in Fig. 1 (N = 6) corresponds at low defor- 
mations 5 — > 1 to a closed shell quantum dot. When analysing the Ml and E2 responses 
of dots that correspond to open shells in the circular limit the results are quite similar to 
those in this Figure, but with the important difference that the u>_ mode is already active 
at B = since Pauli blocking is not effective in this case. 



III. THE LSDA APPROACH 
A. The method 

To analyze the role of electron-electron interactions in a microscopic formalism we resort 
to the LSDA version of density functional theory. This theory provides a practical way to 
include electronic exchange and correlation effects by relying on exact calculations for the 
uniform electron gas. At different levels, density functional theory has been applied by many 
authors to describe the ground state p5|-p^] and excitations of quantum dots ||7|-P,pO|. The 
most refined version is the so-called current-density functional, first applied to quantum dots 



by Ferconi and Vignale [25|, that was recently used to describe the edge reconstruction in 



these systems for large magnetic fields pT|. Nevertheless, current terms are known to be 
rather small for moderate magnetic fields and will be neglected in this work in which we shall 
resort to the approach of Ref. ||, developed for the treatment of noncircular nanostructures. 
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The description of excitations in deformed nanostructures constitutes a highly nontrivial 
task, mainly because of the lack of symmetry not allowing the analytic integration of angular 
variables as in circular systems. For this reason we shall use the time- dependent LSDA to 
obtain the time evolution, following an initial perturbation, of relevant expectation values. 
The set of single-particle orbitals {<pj(r)} evolves in time as 

d 

^V^(r, t) = h v [p, m) cp iri (r, t) , (16) 

where rj =1,1 is the spin index, and total density and magnetization are given in terms of 
the spin densities p n {v) = Yh I^MI 2 ) by p = pj + p\ and m = p-f — p|, respectively. The 
Hamiltonian h v in equation fll6|) contains, besides of kinetic energy, the confining v^ con ^(r), 
Hartree v^ H \r) = J dr'p(r')/|r — r'| and exchange- correlation v^ xc \r) = ^-£ xc (p,m) po- 
tentials. The exchange-correlation energy density S xc (p,m) has been described as in Refs. 
|7|-p|,|30|. This technique was already used by us in Ref. |0| to describe orbital (a>_) excita- 



tions in quantum dots at zero magnetic field. 

An initial perturbation of the orbitals <p'(r) = "P<p(r) models the interaction with the 
external field. The unitary operator V is given in terms of a displacement field u(r) as 

V = exp [m(r) • p] . (17) 

The appropriate fields for the Ml and E2 channels are u M1 = Xre g and u E2 = W(xy), 
where A is a small parameter that guarantees the linear regime. Actually, these opera- 
tors correspond to a rigid rotation (Ml) and a quadrupole distortion (E2) of the electronic 
orbitals. After the initial distortion we keep track of (J2i^z){t) or (Y,i x il/i)(t) an d a subse- 
quent frequency analysis of the signal provides the absorption energies and their associated 
strengths M. 



B. LSDA results and discussion 

Figure 3 summarizes the Ml and E2 spectra of the N = 6 electron dot within time- 
dependent LSDA. Quite remarkably many features of the NIDHO model are also present 
in the LSDA spectra. In fact we can identify the four modes, labeled in Fig. 3 by analogy 
with the NIDHO results, although in general they lie at different energies because of the 
interactions. 

Focusing on the Ml channel, we notice for 5 = 0.9 the emergence of the low energy mode 
oj_ and its crossing with the 0*22 mode. This crossing takes place at B pa 1.5 T for both 
5 = 0.9 and 5 = 0.7 and at B pa 1 T for 5 = 0.4. The dominance of U22 at large -B's is also 
in very good agreement with the NIDHO prediction. Quite remarkably, for 5 = 0.7 the u- 
mode is switched on even at B = as a result of the interactions. At 5 = 0.4 the closeness 
of 

(J22 and uj_ at low magnetic fields predicted in the non-interacting case is nicely seen in 
Fig. 3. It is also worth to mention the strong quenching of the u\\ mode within LSDA. 

The E2 results also reflect to a great extent the systematics predicted by the NIDHO 
model. In general the high energy modes u + and u>n are more important in the E2 spectra 
than in the Ml ones, especially at large magnetic fields. In some cases the LSDA strengths 
show some small fragmentations that we attribute to ph effects (Landau damping). For 
instance, this is quite clear in the Wn mode at 3 and 4 teslas. 
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IV. SUMMARY AND CONCLUSIONS 



The Ml and E2 channels of elliptic quantum dots are characterized by four modes whose 
energies and strengths vary with magnetic field and deformation. Their main characteristics 
are already seen in a non-interacting model, although interactions shift the mode energies and 
introduce several minor differences. The main trends seen in Figs. 1-3 may be summarized 
as follows: 

a) Three modes Wn, u + and uj_ have a positive dispersion relation with B while the other 
U22 exhibits a negative one. At low magnetic fields and deformations wn, uj + and U22 
are very close, while uj- lies at a much lower energy. As the deformation is increased 
the three upper modes separate in energy (at B & 0), the lower of them U22 coming 
close to c<j_. 

b) At low B's the Ml strength of closed shell dots is exhausted by uo + while, when 
increasing the deformation u_ is also active. This mode is also switched on when 
increasing the magnetic field for a fixed deformation as a consequence of the Pauli 
blocking mechanism. Generally, a crossing between u_ and uj 2 2 occurs for intermediate 
magnetic fields. At large magnetic fields the low energy cj 2 2 mode eventually takes all 
the strength. 

c) The E2 spectra are similar to the Ml ones, although an important amount of strength 
shifts to the high energy modes u+, ou- and wu as compared to the magnetic dipole 
channel. Contrary to the Ml case the E2 absorption does not show a clear dominant 
peak at high 5's. 
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FIGURES 



FIG. 1. Results within the NIDHO model for a N = 6 electron dot with r s = 1.51. Left panels 
display the mode energies in H* while central and right ones show the Ml and E2 absorption 
cross sections ujS(uj) as a function of the magnetic field. Each row corresponds to a different dot 
deformation. 

FIG. 2. The lines show the boundary between regions for u>- mode in the NIDHO model and 
different dot sizes. Left and right of each curve correspond to U- being active or non-active (Pauli 
blocked) respectively. 

FIG. 3. Ml an E2 absorption intensities within LSDA for different magnetic fields and defor- 
mations. The results correspond to a dot with N = 6 electrons and r s = 1.51. 
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